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ABSTRACT 

The rich SMC star cluster NGC 419 has recently been found to present both a 
broad main sequence turn-off and a dual red clump of giants, in the sharp colour- 
magnitude diagrams (CMD) derived from the High Resolution Channel of the Ad- 
vanced Camera for Surveys on board the Hubble Space Telescope. In this work, we 
apply to the NGC 419 data the classical method of star formation history (SFH) 
recovery via CMD reconstruction, deriving for the first time this function for a star 
cluster with multiple turn-offs. The values for the cluster metallicity, reddening, dis- 
tance and binary fraction, were varied within the limits allowed by present observa- 
tions. The global best-fitting solution is an excellent fit to the data, reproducing all 
the CMD features with striking accuracy. The corresponding star formation rate is 
provided together with estimates of its random and systematic errors. Star formation 
is found to last for at least 700 Myr, and to have a marked peak at the middle of this 
interval, for an age of 1.5 Gyr. Our findings argue in favour of multiple star forma- 
tion episodes (or continued star formation) being at the origin of the multiple main 
sequence turn-offs in Magellanic Cloud clusters with ages around 1 Gyr. It remains to 
be tested whether alternative hypotheses, such as a main sequence spread caused by 
rotation, could produce similarly good fits to the data. 
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1 INTRODUCTION 

Back in the nineties, the Hubble Space Telescope Wide Field 
Planetary Camera 2 (HST/WFPC2) opened a new era in 
the study of s tellar population s in th e Magellanic Clouds. 
Starting from iGallagher et alj (1 19961 ). many authors were 
able to derive the detailed star formation history of LMC 
and SMC fields, via the analysis of deep CMDs reaching 
well below the oldest main sequence turn-offs (MSTO). It 
became possible as well the accurate measurement of the 
properties and structure of many populous star clusters (e.g . 
Mighell et al.lll996l;lElson et alj 19981; iRich et alj|200d . l200lj ; 



Gouliermis et alj|2004 iKerber et al.ll2007f ) 



The higher efficiency and larger area of the Advanced 
Camera for Surveys (ACS) further improved the situation. 
One of the latest achievements in the field was the con- 
clusive evidence, based on ACS data, that many star clus- 
ters in the LMC with ages typically larger than 1 Gyr, 
presen t double or multiple main sequence turn-offs (MM- 
STO; iMackev fc Brobv Nielsenl 120071 ; iMackev et all |200S|; 



iMilone et ai1l2009r ). This effect was, in precedence, already 
indicated by groun d-based deep CMDs (|Bertelli et al.ll2003l ; 
iBaume et al.ll2007l ) , which however were far from presenting 
the fine details of the ACS ones. 

While there is firm observational basis for the presence 
of MMSTOs in star clusters, their interpretation is far from 
being settled. Once it has been demonstrated that the pres- 
ence of bin aries cannot mimick the detailed shape of MM- 
STOs (see IMackev et al.ll20Qg| ; iGoudfrooii et al.[|2009l ) the 
most obvious interpretation is that they are the signature 
of stellar populations spanning several hundreds of Myr in 
age. This interpretation however poses a major challenge to 
the understanding of star formation and dynamics in star 
clusters, since it is not obvious how objects with relatively 
shallow potential wells could retai n their gas and c o ntinue 
forming stars for so long a time. iBekki fc Mackevl (|2009T l 
propose a mechanism to explain the bluest MSTOs as being 
due to a second event of star formation driven by the colli- 
sion with a giant molecular cloud; this hypothesis however 
does only explain double MSTOs and not the continuous 
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MMSTO structures observed by iGoudfrooii et all (|2009l) in 

NG C 1846. 

iBastian fc de Minkl (|2009r i have recently advanced a 
completely different explanation for MMSTOs: the phe- 
nomenon could be caused by the effects of rotation during 
the main sequence evolution of stars. For a coeval popula- 
tion, the most rapid rotators would evolve to lower T e g and 
generate the redder turn-offs. The hypothesis is certainly in- 
teresting, although it requires an ad hoc decrease of rotation 
velocities with the stellar mass, and relatively high rotation 
velocities - nearly half the critical break-up rotation rate 
- in order to explain the observed features. Whether these 
rotation velocities can be typical of A and F stars in the 
Magellanic Clouds, and cause the large spread in T e a that 
is apparently observed, is yet to be clarified by new obser- 
vations. 

Whether the MMSTO phenomenum can be explained 
by stellar populations of single or multiple ages, bears very 
much in the interpretation of CMDs of nearby galaxies in 
general. In deriving SFH of galaxies, one assumes that their 
CMDs are made by the superposition of populations of dif- 
ferent ages, each one presenting a narrow MSTO. Were the 
MSTOs of single-burst populations intrinsically broad, most 
works of SFH-recovery in nearby galaxies would have to be 
revised to some extent. 

Given the above situation, we decided to test if the 
extended-SFH hypothesis could really explain the MMSTOs 
in ~ 1-Gyr old star clusters in the Magellanic Clouds. As a 
first target, we choose the rich SMC star cluster NGC 419. 
Two aspects make this cluster an ideal ta rget for our goals . 
First, it does not only present a MMSTO rtGlatt et alj|2008h 
but a lso has a dual clump of red giants ( Girardi et alj 
l2009al ). which is the signature of stars close to the tran- 
sition between those that form an electron-degenerate core 
after H-exhaustion, and those that do not. This additional 
fine structure of the CMD provides strong constraints to 
the cluster mean age a nd to the efficiency o f convective core 
overshooting in stars (|Girardi et al.|[2009al ). probably help- 
ing to limit the family of stellar models that can be fit in 
the process of SFH-recovery. Second, the central region of 
NGC 419 counts with extremely accurate photometry pro- 
vided by the High Resolution Channel (HRC) of ACS. As a 
bonus, the field contamination is almost negligible for this 
cluster. 

This paper is structured as follows: Sect. [5] briefly de- 
scribes the SFH-recovery method and its application to the 
NGC 419 data. Sect. |3] presents the results for NGC 419 and 
a discussion about the random and systematics errors in the 
recovered SFH, due to the statistical fluctuations and uncer- 
tainties in reddening, distance, metallicity, binary fraction, 
and field contamination. Sect. U draws the final conclusions. 



2 PREPARING NGC 419 FOR 
SFH-RECOVERY 

2.1 Data and photometry 

The dataset used in this paper has been retrieved from the 
HST archive (GO-10396, PI: J.S. Gallagher) and consists 
of a 740 arcsec 2 area centered on NGC 419, observed with 
the ACS/HRC in the filters F555W and F814W. The same 
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Figure 1. The CMD for NGC 419 as derived from the HRC 
data centered on the cluster, after the re-reduction described in 
this paper and without applying any quality cut to the photom- 
etry. The la error bars, as derived from artificial star tests (see 
Sect. [2721 . are drawn at the left. 



imag es ha ve already been analy s ed by iGlatt et al.l (I20081 . 



ages na ve already b een analysed by mat t et al.l {\ZU 
120091 ) and I Girardi et ail (|2009aE . iGirardi et all (|2009bl 7 its 



figure 1) provide a false-color version of these HRC images, 
together with a comparison with the much wider ACS Wide 
Field Channel (WFC) images of the same cluster. There, 
one may notice that the HRC field samples just the very 
inner part of the cluster, and that the distribution of stars 
over its area is quite uniform. 

The origina l data has been re-reduced as described 
in IGirardi et al.l (|2009al ) . Essentially, we have accessed the 
archival images already processed a nd calibrated u si ng th e 
standard procedures mentioned in ISirianni et al.l (|2005h . 
Aperture and Point Spread Function (P SF) photometry 
with the DAOPHOT package (|Stetsonl 19871 ) were performed 
within a 2 pix radius, and aperture corrections were applied. 
The resulting CM Ds are very similar to those described in 
Idatt et all (|200Sl ). The PSF photometry was then preferred 
and used for all subsequent applications. 

Fig. Q] shows the ACS/HRC data for NGC 419, in the 
F814W vs. F555W-F814W CMD. This plot will be used 
as a reference to our analysis. 



The SFH of NGC 419 3 




0.9 



0.8 



0.0 0.5 1.0 

F555W-F814W 



1-0.4 
™0.3 



Figure 2. Completeness map, derived from the complete set of 
ASTs. 



2.2 Assessing photometric errors and 
completeness 

In order to characterize the errors in the photometry and 
the completeness of the sample, we have performed a series 
of artificial star tests (AST) on the reduced images (see e.g. 
iGallart et aT]|l999l ; lHarris fc Zaritsky||200 J ) . The procedure 
consists of adding stars of known magnitude and colour at 
random places in the F555W and F814W images, and redo- 
ing the photometry exactly in the same way as before. The 
artificial stars are then searched in the photometric cata- 
logues, and when recovered the changes in their magnitudes 
are stored for subsequent use. 

In order to avoid the introduction of additional crowd- 
ing in the images, artificial stars are positioned at distances 
much higher than their PSF width. We found that the PSF 
radius in on our ACS/HRC images is of < 7 pix, whereas 
our fitting radius is of 2 pix. So, our AST are distributed 
on a grid spaced by 20 pix, which is each time randomly 
displaced over the image. 

A total of 3.4 x 10 9 ASTs were performed, covering in 
an almost uniform way the CMD area of the observed stars 
as well as the area for which we build the "partial models" to 
be used in the SFH analysis (see Sect. l2.3l below). Then, the 
ratio between recovered and input stars gives origin to the 
completeness map of Fig. [2] Note that the 90 % completeness 
limit is located at F814W~ 21.7, which is well below the 
position of the MMSTOs in NGC 419. 

Another important aspect is that the stellar density is 
quite constant over the HRC images, as well as the com- 
pleteness. For instance, at magnitudes 22.15, from the image 
center to the borders the completeness changes from 0.87 to 
0.90 for the F555W filter, and from 0.78 to 0.85 for F814W. 
Since this difference is very small and the bottom part of the 
CMD will not be used in our analysis, we do not apply any 
position-dependent completeness to the CMDs, but simply 
the average values derived from all the ASTs. 
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Figure 3. Map of photometric errors as a function of input 
F555W and F814W, as derived from the ASTs. The errors arc 
defined as the difference between the recovered and input magni- 
tudes. 



Figure [3] presents the differences between the recovered 
and input magnitudes of ASTs, as a function of their in- 
put magnitudes. These differences give a good handle of the 
photometric errors effectively present in the data. The error 
distributions are slightly asymmetric because of crowding. 



2.3 The partial models 

The basic process of SFH-recovery is the decomposition of 
an observed CMD as the sum of several "partial models", 
which represent stellar populations in very limited intervals 
of age and metallicity. In our case we will assume a con- 
stant age-metallicity relation (AMR) i.e., a single value for 
the metallicity for all ages, since so far there are no evi- 
dences for signific ant spreads in metallici ty in star clusters 
like NGC 419 (e.g. lMucciareiu" et al.|[2008l ). Hence, each par- 
tial model is defined as a stellar aggregate with constant star 
formation over an age interval of width Alogt = 0.05 dex. 
This is a fine age resolution for a work of SFH-recovery; 
suffice it to recall that the age bins adopted in the SFH- 



in log t (see e. 


y Gallart et al.lll999l; lHarris & Zaritskvll200ll; 


iDolohin et al. 


20031. for some examples). The age interval 



covered by our partial models goes from log(f/yr) = 8.9 
to 9.4, which is much wider than the interval suggested by 
the position of NGC 419 MMSTOs. So, for each set of pa- 
rameters, we have a total of 10 partial models, completely 
encompassing the age interval of interest. 

The other parameters defining a set of partial mod- 
els are the V-band extinction Av, the distance modulus 
(m — M)o, the metallicity [Fe/H], and the binary fraction /b; 
these will be varied in our analysis as described below. Other 
parameters describe the area of the CMD to be sampled, and 
its resolution. With basis on our Figs. [1] and [2] we decide to 
limit the CMD area to be analysed as: F555W - F814W be- 
tween -0.4 and 1.6, and F814W between 16.32 and 22.0. 
Within these limits, stars are not saturated, and complete- 
ness is above 80 %. The CMD resolution is set to be 0.02 mag 
both in colour and in magnitude. 

The partial models are simulated with t he aid of the 
stellar population synthesis tool TRILEGAL (|Girardi et alj 
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Figure 4. Hess diagram for a single partial model, before (left 
panel) and after (right panel) applying the results of ASTs. The 
partial model has a mean age log(t/yr) = 9.125, metallicity 
[Fe/H] = —0.95, and binary fraction /b = 0.18. The colour scale 
indicates the density of stars. 

120051) in its version 1.3, which stands on the iMarigo et all 
(|2008f ) isochrones, transformed t o the A CS/HRC Vegamag 
photometry using iGirardi et al.l (|2008l ) bolometric correc- 
tions a nd extinct i on co efficients. The IMF was assumed to 
be the IChabrierl (|200ll ) lognormal one. The output cata- 
logues from TRILEGAL arc first degraded applying the re- 
sults from the ASTs - which results in the blurring and 
depopulation of the bottom part of the CMDs - and then 
converted into Hess diagrams. Fig. [4] illustrates one of such 
diagrams for a partial model before and after the results of 
ASTs have been applied. 



2.4 The range of cluster parameters 

In our analysis, we have tried to cover the complete possible 
range of parameters for NGC 419, in particular: 



needed. In this work, we have decided to explore the whole 
-0.5 < [Fe/H] < -1.1 interval. 

Binary fraction fa: lElson et al.l (|l99ct ). using the shape of 
the main sequence in a deep CMD (obtained with HST) for 
the LMC cluster NGC 1818, determined a binary fraction 
from ~ 0.20 to ~ 0.35, from the center to the outer parts of 
the cluster. These binary fractions refer only to the binaries 
with high primary /secondary mass ratio, say of above 0.7, 
because these are the only ones which separate clearly from 
the single-sta r main sequence in CM Ds. In our previous work 
on NGC 419 jGirardi et al.ll2009al ) , we have already noticed 
that the HRC CMD was indicating a binary fraction of the 
order of 0.2, consistent with the central region of NGC 1818. 
In the following, we will adopt a reference value of fa ~ 0.18; 
notice however that we have explored fa values from 0.10 to 
0.28. The distribution of primary /secondary mass ratios is 
assumed to be uniform over the interval from 0.7 to 1.0 

Compared to the other parameters, the binary fraction is 
the less important in defining the quality of a CMD fit, 
as will be illustrated below. A clear advantage of explor- 
ing ranges for the other parameters - instead of taking sin- 
gle values from the literature - is that in this way we can 
partially compensate for the possible errors and offsets in 
the theoretical models: Errors in the zeropoints of bolomet- 
ric corrections and on the helium content of stellar models 
could be easily compensated by a change in apparent dis- 
tance modulus (m — M)o + Av- Errors in the theoretical T e g- 
color relations can also be compensated by small changes in 
both Ay and [Fe/H]. Therefore, the final best-fitting values 
for these quantities are to be associated to the stellar models 
which we are using. Different sets of stellar models are likely 
to produce slightly different best-fitting values. 



Distance modulus (m— M)o'- Our reference value for 
the SMC distance modulus is (m-M)o = 18.89 ± 0.14 

|Harries et al.1 120031 ). which is a valu e very simi lar to 

many other recent determinatio ns^. g. ICrowl et alj l200ll : 
lLah et al.ll2005l ; iNorth et al.ll2009T ) . In practice, however, we 
explored (m — Af)o values between 18.75 and 18.97. 



Extinction Ay: The lSchlegel et all (|l998h extinction maps 
provide Eb-v — 0.10 ± 0.01 for an area of radius equal to 5' 
around NGC 419, which transla tes into Ay = . 32 ±0 .03. A 
similar search in the maps from IZaritskv etHI l|2002l) . with 
a 12' radius, produces Ay — 0.27 ± 0.28. In practice, we 
explored Ay values between 0.12 and 0.38. 



Metallicity 

iDurand et al 



l|Kav ser et al 



Fe/H] : Works based on isochrone fitting (e.g. 



1984) as well as Ca n triplet observations 
2009T ) have suggested [Fe/H] values similar to 



— 0.7±0.3 for NGC 419. Interestingly, the metal abundances 
of other SMC clusters with similar ages are uncertain as 
well. Suffices it to mention t he rec ent AMR for the SMC 
as derived from IParisi et ail (|2009l . their figure 14), which 
shows SMC clusters with ages smaller than 2 Gyr covering 
the complete [Fe/H] interval between —0.5 and —1.1. The 
situation is far from clear and high resolution spectroscopy 
of a handful of giants in 1-Gyr old SMC clusters is definitely 



3 RECOVERING THE STAR FORMATION 
HISTORY 

3.1 Method 

To recover the SFH of NGC 419, we use an adapted ver- 
sion of the pipeline that is being built to analyse data from 
the VISTA surv ey of the Magellanic Clouds (VMC; see 
ICioni et al.ll2008T ) . The meth od is fully desc r ibed, and tested 
on VMC simulated data by iKerber et al.1 (|2009l V In brief, 
after building the Hess dia gram for the data and partia l 
models, the StarFISH code (|Harris fc Zaritsky|[200ll , 120041 ) 
is used to find the linear combination of pa rtial mod e ls tha t 
minimizes a \ -like statistic as defined in iDolphinl (|2002h . 
The solution is characterized by the minimum \ 2 , Xmim an d 
by a set of partial model coefficients corresponding to the 
several age bins. These latter translate directly in the star 
formation rate as function of age, SFR(t). 

The SFH-recovery is repeated for each of the (m — M)o, 
Ay, [Fe/H], and fa values. In order to limit the space of 
parameters to be covered, the procedure is essentially the 
following: for a given value of [Fe/H] and /b, we cover a sig- 
nificant region of the (m — M)o versus Ay plane, performing 
SFH-recovery for each point in a grid, and then building a 
map of the Xmin f° r the solutions. Examples of these maps 
are presented in Fig. [5] The maps are extended so that the 
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minimum Xmin f° r a given value of [Fe/H] can be clearly 
identified, as well as the regions in which Xmin increases by 
a factor of about 1.5. The typical resolution of such maps is 
of 0.01 in (m-M)o and 0.005 in Ay. 

The procedure is done for several values of [Fe/H]. Xmin 
maps for different [Fe/H] present relative minima in slightly 
different positions, so that the final resulting maps do not 
cover exactly the same Ay and (m — M)o intervals. In any 
case, the relative minima are well delimitated in all cases, 
as illustrated in Fig. [8] The [Fe/H] spacing between these 
maps is of A [Fe/H] = 0.07 dex. 

Among these many SFH-recovery experiments, the 
most interesting ones are obviously those with the relative 
minima and the absolute minimum Xmin values. Fig.[6]shows 
the x 2 and residuals in the Hess diagram for the best fitting 
solution, which is for (m — M)o = 18.83, Ay — 0.345 and 
[Fe/H] = -0.88 and Xmin = 0.933. It is already evident, 
in this plot, that the fitted solution is an excellent repre- 
sentation of the observed data, with residuals uniformly dis- 
tributed over the CMD. The only point of the CMD in which 
there seems to be some increased residual is the bottom part 
of the main sequence, where according to Figs. [5] and the 
completeness is significantly smaller than 1 (although by just 
10 to 20 %), and photometric errors are higher. 

Figure [7] shows two examples of recovered SFR(t) for 
best-fitting solutions obtained for different values of [Fe/H]. 
They are qualitatively similar, with a clear presence of 
stars spanning ages from ~ 1.2 Gyr (log(t/yr) = 9.08) to 
~ 1.9 Gyr (log(t/yr) = 9.28). 



3.2 Evaluating the errors 

To evaluate the errors for all involved parameters the first 
step was to find the correspondence between the Xmin value 
for each model and its significance (or confidence) level, a. 
This correspondence was estimated simulating 100 synthetic 
CMDs generated with a number of stars equal to the ob- 
served CMD, using the best-fitting SFR(t) and its param- 
eters as the input for the simulations. So, after recovering 
the SFH for this sample of synthetic CMDs, it was possible 
to build the Xmin distribution and to establish the relation 
between the Xmin difference above the minimum and a. 

Figure [8] shows the significance levels for all the solu- 
tions depicted in Fig. [5] As one sees, for all [Fe/H] between 
—0.74 and —0.95 there are ample areas of the Ay versus 
(m— M)o diagrams with solutions within the 95 % signifi- 
cance level. The best solutions, with Xmin close to 0.95 and 
significance levels of about 30 %, are limited to the central 
regions of the [Fe/H] = —0.81, —0.88, and —0.95 diagrams. 

These maps can be used to estimate the mean values of 
the parameters that characterize the best-fitting solutions, 
and their errors. These are determined as the mean and 
standard deviation, weighted by 1 — a, inside the regions in 
which a < 0.68. The results are: 



(m-M)o 
Ay 
[Fe/H] 



18.84 ± 0.04 
0.33 ±0.05 
-0.86 ±0.09 




These numbers give an idea of the region of parameter space 
actually covered by the best-fitting solutions. 

Fig. [9] presents the SFR(t) for the overall best fitting 



9.2 
log(r/yr) 

Figure 9. The overall best-fitting solution for the SFR(f) of 
NGC 419. The random (error bars) and systematic errors (solid 
thin lines) are also shown. 



solution, obtained at [Fe/H] = —0.88, and corresponding to 
the Hess diagrams of Fig. [6] Overplotted on the SFR(t), we 
have the random errors due to the statistical fluctuations 
(associated to the number of stars); they were determined 
from the r.m.s. dispersion in the recovered SFR(i) for the 
100 synthetic CMDs. On the other hand, the systematic er- 
rors in the SFR(t) due to the uncertainties in (m — M)o, 
Ay and [Fe/H], were determined using the minimum and 
maximum SFR(t) values for all models inside the 68 % sig- 
nificance level for all metallicities. 



3.3 The role of the binary fraction 

As anticipated previously, we have verified that the binary 
fraction /b plays a relatively minor role in determining the 
best-fitting solution and the SFR(£). Many simulations at 
varying /b were tried at several stages of our analysis, al- 
ways supporting this conclusion, and always pointing to a 
binary fraction close to 0.2. Fig.[l0]aims to illustrate the sort 
of results one gets for different /b. To build the figure, we 
first define a large set of SFH-solutions for which the CMD 
fitting was considered to be good: namely the entire 68 % 
significance level area for the metallicity [Fe/H] — —0.81, in 
Fig- [HI For this region of the Ay versus (m — M)q plane, the 
SFH-recovery tests are repeated for all /b values between 
0.10 and 0.28, at steps of 0.02. Then, we plot the value of 
mean Xmin an d its r.m.s. dispersion as a function of /b, as 
shown in Fig. 1101 The minimum of mean Xmin f° r /b = 0.18 
is evident. However, the dispersion of Xmin values is also 
significant, and indicates that a relatively wide range of /b, 
from say 0.12 to 0.24, would also produce acceptable results, 
if compared to the fh = 0.18 case. 

It is also worth noticing that for /b = 0.18 the 68 % 
significance level (a = 0.68) is about 0.038 above the mini- 
mum Xmini as determined in the previous section. All solu- 
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Figure 6. In the Hess diagram, we show the data (left), the solution found by StarFISH (middle), the data-model difference and the x 2 
map. The model is for [Fe/H] = -0.88, (m-M) = 18.83, A v = 0.345, f h = 0.18. 

tions presented in Fig. [10] are inside this limit. So, even if to the random errors in Fig.[9j but they also do it in a non- 
there is a apparent minimum Xmin around /b = 0.18, it is systematic way: the mean age and shape of the SFR(t) do 
not statistically significant. not show any significant trend with /b. 

Finally, Fig. Illl shows the recovered SFR(t) for four dif- 
ferent fh values over the /b = 0.12 to 0.28 interval. Not only 
the SFR(t) is seen to vary by amounts that are comparable 
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Figure 7. The SFR(t) corresponding to the best fitting solution (blue histogram), for metallicitics [Fc/H] = —0.95 (left) and —0.81 
(right). The error bars indicate the random errors, whereas the green histograms indicate the systematic errors (see Sect. IQfor more 
details). 
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Figure 8. Significance level distributions for several [Fe/H] values. 
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Figure 10. The average of the minimum Xmin vers us the binary 
fraction /b, for a series of models of metallicity [Fc/H] = —0.81. 
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Figure 11. The SFR(t) for a series of SFH solutions with 
[Fe/H] = -0.81, (m-M) = 18.86, A v = 0.315, and for dif- 
ferent binary fractions: /j, = 0.12 (blue dashed line), 0.18 (red 
dashed), 0.24 (blue solid), and 0.28 (red solid). 



3.4 The role of the field contamination 

In this work, we have analysed the ACS/HRC photometry 
NGC 419 without taking into account that a fraction of the 
observed stars is due to the SMC field. Indeed, when first 
notici ng the unusual str ucture of the red clump in these 
data, Idatt et al.l (|2008l ) suggested the fainter red clump 
could be due to the SMC field. The counter-argument by 
iGirardi et all (|2009al l . however, has dramatically redimen- 
sioned the possible role of the SMC field contamination in 
this case: after measuring the density of red clump stars in 
the external parts of the ACS/WFC images of the same 
cluster, at radii larger than 75" and for a total area of 
2.47 x 10 4 arcsec 2 , they find that just 4.5 red clump stars 
from the SMC field are expected to be found inside the 
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Figure 12. The SFR(i) and AMR obtained when the [Fe/H] 
is not constrained to a single value, but is let to vary between 
4 different values in the -0.95 ^ [Fe/H] ^ -0.80 interval. The 
upper panel shows the [Fe/H]-added SFR(t), whereas the bottom 
panel show the SFR-averaged [Fe/H](t). 



740 arcsec 2 area of the HRC images, whereas the total ob- 
served number is of 388 (47 of them are in the secondary 
clump). These numbers set the probability that the sec- 
ondary clump i n NGC 419 is made of SMC field stars to 
less than 1(T 9 (|Girardi et al.ll2009bl ). Moreover, they indi- 
cate that the field contamination in the red clump area of 
the CMD is of just ~ 1 %. 

In the context of the SFH-recovery work, some addi- 
tional numbers worth of consideration are the following. The 
number of field main sequence stars in the magnitude inter- 
val 20 < F814W < 22, for the same 2.47 x 10 4 arcsec 2 area 
in the outskirts of the WFC images, is 2640, whereas the 
observed number in the 740 arcsec 2 of HRC is 2395. This 
magnitude interval is just barely affected by incompleteness, 
and contains most of the observed stars. Therefore, we can 
estimate that just ~ 3.3 % of the stars used to study the 
SFH of NGC 419 are due to the SMC field. This fraction 
is similar to the Poisson noise in the total number of HRC 
stars, indicating that our SFH results cannot be affected in 
a significant way. 



3.5 The age— metallicity relation 

With star formation lasting for 700 Myr, one may wonder 
whether NGC 419 would not have enriched itself with metals 
ejected from its oldest stellar populations. Although we be- 
lieve that the present photometric data is not enough to clar- 
ify this question, we performed an additional SFH-recovery 
exercise, in which the cluster [Fe/H] is not fixed, but can take 
four different values at each age: [Fe/H] = —0.95, —0.90, 
—0.85, —0.80. In practice, this is obtained by running the 
SFH-recovery method with a larger library of partial models, 
with 40 components (10 age bins times 4 metallicity bins). 
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This larger set of partial models provides more free- 
dom to the StarFISH algorithm, which converges to an 
absolute best-fitting solution with a Xmin 01 just 0.84, for 
A v = 0.33 and (m-M) = 18.83. The SFH results are 
illustrated in Fig. 1121 the upper panel shows the SFR(t) 
obtained by adding the SFR values obtained for all metal- 
licities at each age bin, whereas the bottom panel shows the 
AMR as obtained by a SFR-weighted average of the four 
different [Fe/H] values. 

This SFR(t) is actually very similar to the one obtained 
form the best-fitting solution with [Fe/H] = -0.88 (Fig. IS]). 
In the age interval with significant star formation, which 
goes from 1.2 to 1.9 Gyr and comprises 5 age bins, the three 
youngest age bins are found to have a metallicity almost 
coincident with the [Fe/H] = —0.88 one, whereas the two 
oldest ones are found to be slightly more metal-rich, with 
—0.80 dex. It is not clear whether this [Fe/H] variation can 
be significant, since it is comparable to the estimated error of 
0.09 dex in the [Fe/H] of the best-fitting solutions (Sect . 15^ . 
We conclude that our analysis does not provide evidence for 
self-enrichment taking place in NGC 419. 



4 CONCLUSIONS 

In this paper, we have derived the SFH of the SMC star 
cluster NGC 419, via the classical method of CMD recon- 
struction. This analysis implicitly assumes that NGC 419 
is formed by a sum of single-burst stellar populations (or 
"partial models"), each one being characterized by a well- 
defined and narrow MSTO. The only effects that blur these 
partial models in the CMDs are the sequences of binaries 
(which effect however is far from dramatic) and the photo- 
metric errors. With this kind of approach, the broad MM- 
STO observed in NGC 419 naturally translate into a con- 
tinued star formation history. We find a SFR(i) lasting for 
700 Myr (from 1.2 to 1.9 Gyr, see Fig. |9j), which is quite a 
long period, probably at the upper limit of all values already 
estimated for star cluster s with MMSTOs (iBertelli et al] 



20031; iBaume et al 



[20071; iMackey fc Brobv Nielsen! 120071; 
Mackev et al.1 120081 ; iMilone et al.l 120091 ; iGoudfrooii et al.1 
20091 ). Our error analysis leaves practically no room for the 



MMSTOs in this cluster being caused by a single episode of 
star formation. 

It is also remarkable that the SFR(£) we derive presents 
a pronounced maximum at the middle of the star forma- 
tion interval, at an age of 1.5 Gyr (Fig. |9}. This is some- 
what unexpected, in the context of the few suggested sce- 
narios for the appeara nce of prolonged star f ormation in 
LMC star clusters (see IGoudfrooii et all 120091 . for a com- 
prehensive di scussion of them). Either the m erging of two 
star clusters l|Mackev fc Brobv Nielsen! 120071 ), or a second 
period of star forma tion driven by the mer ging with a gi- 
ant molecular cloud (|Bekki fc Mackevi[2009h . would lead to 
strongly bimodal distributions of cluster ages, which we do 
not find in NGC 419. Only in the case of the best model 
with [Fe/H] = —0.95, there is an indication of two differ- 
ent peaks in the SFR(f) (see Fig. [7J, which however are 
neither separated nor followed by periods of null SFR(f). 
Also the trapping of field stars by t he forming star cluster 
|Pflamm-Altenburg fc Kroupall2007h would lead to different 
results, with the major peak of star formation being found at 



the youngest ages. O n the other hand, our fi ndings seem to 
be more in line with IGoudfrooii et al.1 HjojJ's conclusions, 
based on the quite continuous distribution of stars across the 
MMSTO region of the LMC cluster NGC 1846. They sug- 
gest a scenario in which the star cluster continues to form 
stars in its center out of the ejecta of stars from previous 
generations. In our case, however, this process would have 
to proceed for a significantly more extended period of time 
than for NGC 1846, and peak - somewhat against the most 
naive expectations - at the middle of the star formation pe- 
riod. Also, this latter scenario might imply some amount of 
self-enrichment in this cluster, whereas our method instead 
is compatible either with a constant metallicity, or with some 
amount of metal dillution. Needless to say, the present ob- 
servational data is not clear enough to provide unambiguous 
indications about the formation scenario of such clusters. 

Another basic result of our analysis is that the hypoth- 
esis of continued SFH, together with current stellar evolu- 
tionary models and a modest fraction of binaries, produces 
excellent fits to the CMD of NGC 419, with Xmin as small 
as 0.93 - or even 0.84 if the [Fe/H] is let to vary during 
the SFH-recovery. We have translated this result into quan- 
titative estimates for the random and systematic errors of 
the derived SFR(t). It is obvious that many combinations of 
cluster parameters produce acceptable solutions, with sig- 
nificance levels smaller than 95 %. However, the really good 
solutions - i.e. those with significance levels better then 68 % 
- cover a narrow region of the parameter space, comprising 
less than 0.1 mag in both (m — M)o and Ay, and about 
0.1 dex in [Fe/H]. 

Despite our success in fitting the CMD of NGC 419 
as a sum of partial models, this success does not tell us 
about the reliability of alternative hypotheses for the origin 
of M MSTOs. In particular, we cannot conclude anything 
about iBastian fc de Minkl (|2009T l 's hypothesis based on the 
colour spread of coeval stars with different rotation rates, 
since its capability of providing a quantitative description 
of the data has not yet been tested. We urge that this test 
should be performed, in order to shed light into this problem. 

However, we call attent ion to another quest i on tha t has 
to be properly addressed: Is lBastian fc de Minkl (|2009T )'s ro- 
tation hypothesis able to produce, in addition to the broad 
turn-off of Magellanic Cloud star clusters, also the dua l red 
clump observed in NGC 419 (see IGirardi et al.ll2009ah ? In 
our models with prolonged star formation, the 700 Myr 
spread in age translates into a ~ 0.26 A'Iq spread in turn-off 
masses, which in turn implies a small spread in the mass 
of H-exhausted cores as stars leave the main sequence. This 
small spread - of the order of just ~ 0.01 Mq - is enough 
to cause the appearance of a dual red clump in t his cluster, 
as thoroughly discussed in IGirardi et al . ( 2009a. b). Such a 
feature is then naturally present in our best-fitting solutions, 
as can be appreciated in Fig. [5] Can rotation do the same, 
producing spreads of a compatible magnitude both in the 
main sequence and in the red clump ? Moreover, can this be 
achieved with a reasonable and simple enough distribution 
of rotation velocities ? At first sight, this seems very unlikely 
to us. 



10 Rubele, Kerber & Girardi 



ACKNOWLEDGMENTS 

We thank the anonymous referee for the useful sugges- 
tions. The data presented in this paper were obtained 
from the Multimission Archive at the Space Telescope Sci- 
ence Institute (MAST). STScI is operated by the Associ- 
ation of Universities for Research in Astronomy, Inc., un- 
der NASA contract NAS5-26555. We thank the support 
from INAF/PRIN07 CRA 1.06.10.03, contract ASI-INAF 
1/016/07/0, and the Brazilian agencies CNPq and FAPESP. 



REFERENCES 

Bastian N., de Mink S. E., 2009, MNRAS, 398, Lll 
Baume G., Carraro G, Costa E., Mendez R. A., Girardi 

L., 2007, MNRAS, 375, 1077 
Bekki K., Mackey A. D., 2009, MNRAS, 394, 124 
Bertelli G, Nasi E., Girardi L., Chiosi C, Zoccali M., Gal- 

lart G, 2003, AJ, 125, 770 
Chabrier G, 2001, ApJ, 554, 1274 

Cioni M., Bekki K., Clementini G, de Blok W. J. G, Emer- 
son e. a., 2008, Publications of the Astronomical Society 
of Australia, 25, 121 

Crowl H. H., Sarajedini A., Piatti A. E., Geisler D., Bica 
E., Claria J. J., Santos Jr. J. F. G, 2001, AJ, 122, 220 

Dolphin A. E., 2002, MNRAS, 332, 91 

Dolphin A. E., Saha A., Skillman E. D, Dohm-Palmer 
R. G, Tolstoy E., Cole A. A., Gallagher J. S., Hoessel 
J. G, Mateo M., 2003, AJ, 126, 187 

Durand D., Hardy E., Melnick J., 1984, ApJ, 283, 552 

Elson R. A. W., Sigurdsson S., Davies M., Hurley J., 
Gilmore G, 1998, MNRAS, 300, 857 

Gallagher J. S., Mould J. R., de Feijter E., Holtzman J., 
Stappers e. a., 1996, ApJ, 466, 732 

Gallart G, Freedman W. L., Aparicio A., Bertelli G, 
Chiosi C, 1999, AJ, 118, 2245 

Girardi L., Dalcanton J., Williams B., de Jong R., Gal- 
lart C, Monelli M., Groenewegen M. A. T., Holtzman 
J. A., Olsen K. A. G, Seth A. C, Weisz D. R., the 
ANGST/ANGRRR Collaboration 2008, PASP, 120, 583 

Girardi L., Groenewegen M. A. T., Hatziminaoglou E., da 
Costa L., 2005, A&A, 436, 895 

Girardi L., Rubele S., Kerber L., 2009a, MNRAS, 394, L74 

Girardi L., Rubele S., Kerber L., 2009b, preprint 
(arXiv:0910.3139) 

Glatt K., Grebel E. K., Gallagher J. S., Nota A., Sabbi 
E., Sirianni M., Clementini G, Da Costa G, Tosi M., 
Harbeck D., Koch A., Kayser A., 2009, AJ, 138, 1403 

Glatt K., Grebel E. K., Sabbi E., Gallagher J. S., Nota A., 
Sirianni M., Clementini G, Tosi M., Harbeck D., Koch 
A., Kayser A., Da Costa G, 2008, AJ, 136, 1703 

Goudfrooij P., Puzia T. H., Kozhurina-Platais V., Chandar 
R., 2009, AJ, 137, 4988 

Gouliermis D, Keller S. G, Kontizas M., Kontizas E., 
Bellas- Velidis I., 2004, A&A, 416, 137 

Harries T. J., Hilditch R. W., Howarth I. D., 2003, MN- 
RAS, 339, 157 

Harris J., Zaritsky D., 2001, ApJS, 136, 25 

Harris J., Zaritsky D., 2004, AJ, 127, 1531 

Kayser A., Grebel E. K., Harbeck D. R., Cole A. A., 
Koch A., Gallagher J. S., da Costa G. S., 2009, The Age- 
Metallicity Relation of the SMC. pp 157 — h 



Kerber L. O., Girardi L., Rubele S., Cioni M.-R., 2009, 

A&A, 499, 697 
Kerber L. O., Santiago B. X., Brocato E., 2007, A&A, 462, 

139 

Lah P., Kiss L. L., Bedding T. R., 2005, MNRAS, 359, L42 
Mackey A. D, Broby Nielsen P., 2007, MNRAS, 379, 151 
Mackey A. D., Broby Nielsen P., Ferguson A. M. N., 

Richardson J. G, 2008, ApJ, 681, L17 
Marigo P., Girardi L., Bressan A., Groenewegen M. A. T., 

Silva L., Granato G. L., 2008, A&A, 482, 883 
Mighell K. J., Rich R. M., Shara M., Fall S. M., 1996, AJ, 

111, 2314 

Milone A. P., Bedin L. R., Piotto G, Anderson J., 2009, 

A&A, 497, 755 
Mucciarelli A., Carretta E., Origlia L., Ferraro F. R., 2008, 

AJ, 136, 375 

North P. L., Gauderon R., Royer F., 2009, in J. T. van 
Loon & J. M. Oliveira ed., IAU Symposium Vol. 256 of 
IAU Symposium, New distance and depth estimates from 
observations of eclipsing binaries in the SMC. pp 57-62 

Parisi M. C, Grocholski A. J., Geisler D., Sarajedini A., 
Claria J. J., 2009, AJ, 138, 517 

Pnamm-Altenburg J., Kroupa P., 2007, MNRAS, 375, 855 

Rich R. M., Shara M., Fall S. M., Zurek D, 2000, AJ, 119, 
197 

Rich R. M., Shara M. M., Zurek D, 2001, AJ, 122, 842 
Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Sirianni M., Jee M. J., Bem'tez N., Blakeslee J. P., Martel 
A. R., Meurer G, Clampin M., De Marchi G, Ford H. C, 
Gilliland R., Hartig G. F., Illingworth G. D., Mack J., 
McCann W. J., 2005, PASP, 117, 1049 

Stetson P. B., 1987, PASP, 99, 191 

Zaritsky D., Harris J., Thompson I. B., Grebel E. K., 
Massey P., 2002, AJ, 123, 855 



